Evidence of Two Distinct Dynamic Critical Exponents in Connection with Vortex Physics 
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The dynamic critical exponent z is determined from numerical simulations for the three-dimensional (3D) 
lattice Coulomb gas (LCG) and the 3D XY models with relaxational dynamics. It is suggested that the dynamics 
is characterized by two distinct dynamic critical indices zq and z related to the divergence of the relaxation time 
r by r oc £ z ° and r oc k~ z , where £ is the correlation length and k the wavevector. The values determined 
are zq ~ 1.5 and z » 1 for the 3D LCG and zq ~ 1.5 and z ~ 2 for the 3D XY model. It is argued that the 
nonlinear IV exponent relates to zo, whereas the usual Hohenberg-Halperin classification relates to z. Possible 
implications for the interpretation of experiments are pointed out. Comparisons with other existing results are 
discussed. 

PACS numbers: 05.70.Jk, 74.40.+k, 75.40.Gb, 75.40.Mg 

A neutral superfluid like 4 He and a superconductor in the 
limit of large London penetration depth can be characterized 
by a complex order parameter. The XY model can be viewed 
as a discretized version of this type of systems in terms of 
the phase of the order parameter [[![]. An interesting feature 
of this class of models is the presence of thermally gener- 
ated topological defects which in three dimensions (3D) take 
the form of vortex loops. The superconducting phase tran- 
sition from the vortex loop point of view separates a low- 
temperature phase consisting of closed vortex loops of finite 
extent from a high-temperature phase where the loops can dis- 
integrate [§-H]. The static thermodynamic properties of the 
XY model allow a dual representation in terms of the vortex 
degrees of freedom. In particular the Villain transformation 
of the XY model gives rise to the dual vortex loop model 
termed the lattice Coulomb gas (LCG) [ph. A crucial point 
in the following is that this duality does not include the dy- 
namic properties which might be quite different. However, 
all these models belong to the same universality class for the 
static phase transition which means that the static critical in- 
dices are the same. 

The universality of the dynamic behavior is weaker and re- 
quires in addition to the static universality also that certain 
global features of the dynamics are the same. In the dynamic 
universality classes defined by Hohenberg and Halperin [^] 
these additional global features are expressed as local conser- 
vation laws. According to this scheme the dynamics of a 3D 
superfluid belongs to model F characterized by the dynamic 
critical exponent z = 1.5 [Q]. A model with purely relax- 
ational dynamics on the other hand belongs to model A with 
z « 2 [H]. In case of a superconductor both model F and 
model A have been proposed as the appropriate dynamical 
class |H|]]. Since relaxational dynamics is related to model A, 
it might appear surprising that the critical dynamic index z for 
the 3D LCG for periodic boundary conditions (PBC) and with 
purely relaxational dynamics in Ref. (verified in Ref. 1 10]) 
was found to be z w 1.5 instead of z » 2. Similarly the 3D 
XY model with fluctuating twist boundary condition (FTBC) 



and relaxational dynamics was found to have z » 1.5 Ql l|]. 
On the other hand, the 3D XY model with PBC and relax- 
ational dynamics has z ~ 2 consistent with model A. This 
implies that the choice of boundary condition affects the value 
of z Jilt], In the present Letter we show that the same type of 
sensitivity applies to the 3D LCG, giving rise to two distinct 
exponents z « 1.5 and z ~ 1. It is proposed that these results 
for the LCG model and the XY model reflect the existence 
of two distinct critical indices corresponding to r oc £ z ° and 
t oc k~ z , where r, £, and k denotes relaxation time, coher- 
ence length and wavevector, respectively. 

The 3D LCG on a cubic lattice is defined by the Hamilto- 
nian [WW^A 
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where = (q a i) represent the vortex line segment variables 
for each lattice site, one for each of the three lattice directions 
e„, q ai £ [0, 1, —1] corresponding to no vortex segment, re- 
spectively to a segment with vorticity 1 or — 1 between neigh- 
boring lattice sites and subject to the constraint that the sum 
of the qai's for the six directed links reaching each lattice site 
is zero. We use PBC and W(r) is the lattice version of the 
Coulomb interaction 11131] 
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where L is the total length of the lattice, the lattice constant 
is 1, and n a = 2sin(A: a /2) where a =1, 2, and 3. The re- 
laxational dynamics is implemented by a Metropolis Monte 
Carlo update where each complete update of the lattice is as- 
sociated with one time unit (for details see Ref. |Q|). The 
voltage across the sample V is proportional to the expan- 
sion rate of the vortex loops V oc JjX)r rxc lr where r 
denotes the site positions on the lattice. The resistance R 
can then be obtained from the voltage fluctuations through 



the Nyquist formula R oc dt(V{t)V(0)}. The scal- 

ing connection, R oc 1/t, leads to the size scaling at T c , 
R oc L~ z °, and the ^-scaling of the resistivity, p oc ^ 1_z ° 
in the critical region [[||. The 3D LCG has a phase transition 
at T c 3.003 with v ps 0.67 Finite-size scaling at T c 

gives zo ~ 1.5 [§Jjjj]- In Fig.l we demonstrate that the size 
converged ^-scaling slightly above T c , using £ oc (T — T c )~ u , 
also gives zq « 1.5. This shows that z$ ps 1.5 is not a finite- 
size effect caused by the boundary, but is a bulk property char- 
acterizing the dynamics. 

The LCG model defined with periodic boundary conditions 
(PBC) corresponds to the XY model defined with the fluctu- 
ating twist boundary condition (FTBC) Jl6|]. The Hamiltonian 
for the latter model is given by 



H(e r ,A) = - 



E 
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where the a summation is over the three nearest neighbors of 
r in the three lattice base directions. The relaxational dynam- 
ics is, as for the LCG model, implemented by a Metropolis 
Monte Carlo update where each complete update of the lattice 
is associated with one time unit (for details see Ref. [|l7[]). We 
use PBC for 9 r which means that A a is the average twist of 
the angle 9(r) across the lattice in the a direction. The fluc- 
tuations in the twist variable A Q (t) are directly related to the 
voltage V a across the sample by V a = — LjjjA a Jis|]. The 
correspondence between LCG with PBC and the XY model 
with FTBC basically hinges on the fact that the models de- 
fined with these boundary conditions allow for voltage fluc- 
tuations across the system. In contrast the usual PBC applied 
to the 3D XY model is equivalent to A = and the volt- 
age fluctuations across the system is then no longer properly 
described by the model [|llj]. The 3D XY model with PBC, 
within the Villain approximation, corresponds to the 3D LCG 
model with PBC described by the modified Hamiltonian [ |l9| | 
[compare Eq.(jl])] 



\ QrtWfa - Yj)q a j + U(2irM a /L 2 ) 
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where 17(0) = YZ=-oo exp[-(0 + 2to) 2 /2T] is the Vil- 
lain function and M a is given by M = |X) r r x < Jr- 
Since this Hamiltonian in Fourier space has the structure 
H K Jlk^o, a < la(^)W(k)q a (~k) + 5wU the difference 
with Eq.(Jl|) is that in the modified Hamiltonian (0) the k = 
mode is suppressed. We focus on the scaling of the vorticity 
correlation function 



G(k,t) = 



k 2 £ 3 



(<?ak(t)?«-k(0)), 
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where k is perpendicular to the a-direction. In particu- 
lar G(k — 0, t) has the scaling form G(k = 0, t) oc 
F(t£,~ z )/£, [pT[]. Figure 2 demonstrates that this scaling is 
well borne out with the value z « 1. The data are from the 
critical region just above T c where the data are size converged 



and £ oc (T — T c )~ v ', which again emphasize, that this is a 
bulk property and not a property that vanishes for L = oo. 
Thus, just as for the 3D XY model in Ref. [p"l|], we obtain 
two distinct values of z pa 1.5 and z » 1, by allowing and 
suppressing the k = mode. We again stress that the scaling 
and the exponents are obtained from data which are indepen- 
dent of the system size. 

As a resolution of this dichotomy we suggest that these 
models are characterized by two distinct indices correspond- 
ing to r oc £ z ° and r oc k~ z . Assuming such a resolution 
would mean that the vorticity correlation function G(k, t) in 
Eq.(|) in general scales as G(k,t) oc F(tk z ,t£~ z ° , 
which precisely at T c for fc m ; n — 2tt/L reduces to 



LG(k min ,t) = F(tk z lin ,tL- Z0 ). 



(6) 



Choosing t = xfc~? n gives F(x,xL z ~ z °) which goes to 
F(x, oo) for large L when z > zq and F(x,Q) for z < zq. 
This means that in the scaling limit of large t (t ^> tq where 
To is a microscopic characteristic time) we will approach a 
scaling F(x) with x = tk^ in for large (small) x when z > zq 
(z < zq). Vice versa we will approach a scaling limit F(x) 
with x = tL~ z ° for large (small) x when zq > z (zq < z). 

We first test the possibility of two distinct indices for the 
3D XY model with FTBC [see Eq.(|)]. The resistance R for 
the 3D XY with FTBC is readily calculated from A a (t) (see 
JT^ ] for details). G(k,t) is obtained by replacing q ar with 
q a (r) = sm(8 r — 8 r+ea ) in Eq.(^]) and multiplying by k 2 



G(k,t) = p(fek(*)ga-k(0)>. 



(7) 



Figure 3(a) and (b) demonstrate, by using the scaling form 
R oc L- Zo f{L{T-T c )- u ) and LG(k,t) = F{tk z min ), that for 
this model zq ps 1.5 whereas z ~ 2, confirming that zq ^ z. 
For Zq < z it should in principle also be possible to observe a 
crossover to F(tL~ z °) for small enough arguments (but still 
with t 3> To), however this limit was not reached in the sim- 
ulations [^T]]. It is interesting to note that the exponent z ps 2 
associated with r oc k* nin for the 3D XY model with FTBC, 
within error bars, has the same value as the exponent zo pa 2 
found for the 3D XY model with PBC [0. This suggests 
that the difference between FTBC and PBC arises from the 
difference in the treatment of the voltage fluctuations across 
the system, or equivalently the k — fluctuations: The k = 
fluctuations for FTBC are associated with zo ~ 1.5. Changing 
to PBC suppresses these fluctuations and the k = fluctua- 
tions for PBC would then corresponds to the fc m i n = 2tt/L- 
fluctuations for FTBC. This argument implies that zo for PBC 
and z for FTBC should be identical, in accordance with our 
numerical data. 

Applying the same reasoning to the 3D LCG with PBC 
would mean that zo ~ 1.5 and z ~ 1. Since in this case 
z < zo, this would mean that LG(k m i n ,t) should scale like 
F(x) with x = k^ iri t for small enough x [compare Eq.(^)]. 
Figure 4(a) shows that this prediction is consistent with the 
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data and that the scaling curve approaches the correct asymp- 
totic form F(x) oc \jx x l z p2| ] in the limit of small x. Figure 
4(b) shows that the data for the largest converged t values are 
instead consistent with the scaling F(x — tL~ Za ), which is 
also in accord with the scaling form given by Eq.(^). From 
this we conclude that the apparent non-uniqueness of the crit- 
ical exponent z for the 3D XY model and the LCG with re- 
laxational dynamics is consistent with the existence of two 
distinct critical indices zq and z for the case when the models 
are defined so as to allow for voltage fluctuations across the 
system. 

We have here proposed the existence of two indices for re- 
laxational dynamics in 3D. On the other hand the correspond- 
ing model with the resistively shunted Josephson junction dy- 
namics (RSJ) has zq — z pl|]22| ]. This model is closely re- 
lated to the XY model and differs in that the dynamics is 
subject to local current conservation, i.e., it has a local con- 
servation law which is not fulfilled by relaxational dynamics. 
In 2D the whole low temperature phase below the Kosterlitz 
Thouless transition is quasicritical and for the 2D XY model 
(with FTBC) Zq(T) = 2irT{T) - 2 > 2 where T is the he- 
licity modulus and this value has been confirmed both for re- 
laxational and RSJ dynamics, whereas z — 2 ^ zq was found 
for the same cases [fit]. Consequently, the existence of two 
indices zq ^ z appears not to be restricted to relaxational dy- 
namics. 

It is interesting to note that the zo value is more universal 
than the z: In 3D z « 1.5 is found for the LCG [|||l0| and 
the XY model (with FTBC) both for relaxational and RSJ dy- 
namics. Likewise in 2D zq(T) = 2irY(T) — 2 is found for 
LCG fl23"| ] (as well as for the 2D Coulomb Gas model with 
Langevin dynamics p4]]) and the XY model (with FTBC) 
both in case of relaxational and RSJ dynamics. On the other 
hand no such universality exists for z: In 3D we have the de- 
creasing sequence z w 2, 1.5, and 1 when going from XY 
relaxational dynamics, XY RSJ dynamics, and LCG relax- 
ational dynamics, respectively. The important point to note 
here is that the choice of dynamics imposed on the vortex lines 
leads to different dynamic characteristics even for the small 
k and long time behavior. One may also observe that the z 
value for the XY model with relaxational dynamics is con- 
sistent with model A value z w 2 in the Hohenberg-Halperin 
classification scheme [|[| which to us suggests that the z de- 
fined from t oc k~ z will in general be consistent with this 
classification scheme. 

A crucial conclusion from the present investigation is that 
which of the two values is appropriate depends on the ex- 
periment performed. For example a very common measure- 
ment for high T c superconductors is the IV characteristics. 
In this experiment the voltage across the sample is measured 
and consequently this measurement relates to zq. The scaling 
prediction for the non-linear IV characteristics gives V oc I a 
where a = (zq + l)/(d — 1) in d dimensions ||]. A compar- 
ison with a = (z + 1) / (d — 1) where z is calculated accord- 
ing to the Hohenberg-Halperin classification scheme is then 
likely to be incorrect if the superconductor is described by 3D 



model A (because zq ^ z for this case), but would be cor- 
rect if it is described by 3D model F (because zq — z for this 
case flTl| , p2| ]). Consequently the possibility of the existence 
of two critical dynamic indices would have to be taken into 
account when analyzing experiments. A criterion for when 
z ^ zq and z = zq remains to be found. 
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FIG. 1. Determination of zo from p oc ^ l zo f° r the 3D LCG 
with PBC for lattice sizes L=10, 12, and 16 (triangles, circles, and 
squares, respectively) and £ oc £ = (T — T c )~ v (see text). As 
seen the data is independent of lattice size. The broken line is a 
least square fit to the linear part of the data and gives 20 ~ 1-4(1) 
consistent with 20 ~ 1.5. The inset shows 2 — 1 from the slope 
(2 — 1 = —slope) for the L — 16-data, obtained by least square fit- 
ting starting from the two data points for the largest £ and then con- 
secutively adding more data points so that the rightmost data point 
in the inset is based on all but one data point in the main figure. The 
dashed lines in the inset give our rough estimate of the error. 




FIG. 2. Determination of 2 from the scaling £G(0, i) oc F(t£ z ) 
for the modified 3D LCG (corresponding to the 3D XY model with 
PBC). The data are for lattice size L = 20 and T = 3.07, 3.10, 3.14, 
and 3.30 (open squares, crosses, filled circles and asterisks, respec- 
tively). A good collapse is obtained for 2 w 1, whereas the inset 
shows that no collapse is obtained for 2 = 1.5 (£ is defined as in 
Fig.l). Consequently z « 1, which is different from z w 1.5 ob- 
tained in Fig.l. 
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FIG. 3. Demonstration that the 3D XY model with FTBC con- 
tains two distinct indices 2 and zo where 2 7^ 20. (a) Scaling of 
resistance, R fa L" z °F(i(T - T c y"), close to and at T c gives 
z w 1.5. (b) Scaling LG(fc min ,i) = F(tk^ in ,tL~ z °) at T c gives 
2 w 2 for large t (L = 10, 12, 16). Note that the scaling with the 
larger of 2 and 20 always dominates in the large t limit. 
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Demonstration that the 3D LCG with PBC at T c contains two distinct indices z and zo within the single scaling function 
, n ,t) = F(tk^ in ,tL~ Z0 ). (a) The data are for L= 6, 8, 10, 12, 16, 20, and 24. The inset shows that LG(k min ,t) does not scale 
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with tfcjnTn 5 for the data shown (the value of LG decreases for increasing L for a fixed tk^^' 5 ). The full drawn part of the data for each size 
in the inset gives a middle section of the data corresponding to neither too large nor too small values of t. These full drawn middle sections of 
the data in the inset are the parts that collapse to a single scaling curve for z m 1, as demonstrated by the main part of the figure. The broken 
curve in the main part is the leading small tfc^ in -dependence of the scaling function given by C/tk^^, where C is a constant. This leading 
term is consistent with the scaling curve obtained from the data (C is used as an adjustable parameter), (b) Data for £=10, 12 and 16 obtained 
for larger t. These larger t data collapse to a single curve with z « 1.4. Note that this is consistent with that the larger of z and zo dominates 
in the large t limit. 
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